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Abstract. The most significant periodicities in the terrestrial 
impact crater record are due to the “hnman-signal”: the bias 
of assigning integer values for the crater ages. This bias seems 
to have eluded the proponents and opponents of real period¬ 
icity in the occurrence of these events, as well as the theorists 
searching for an extraterrestrial explanation for such periodic¬ 
ity. The “human-signal” should be seriously considered by sci¬ 
entists in astronomy, geology and paleontology when searching 
for a connection between terrestrial major comet or asteroid 
impacts and mass extinctions of species. 

Key words: Solar system: Earth - Comets, Techniques: sta¬ 
tistical, Galaxy: solar neighbourhood 


1. Introduction 

An outstanding series of papers appeared in 1984 when a 
28.4 Myr cycle was detected in the terrestrial impact crater 
record (Alvarez & Muller 1984, Davis et al. 1984, Whitmire & 
Jackson 1984). This value was close to the 26 Myr cycle dis¬ 
covered in the geological record of major mass extinctions of 
species (Raup & Sepkoski 1984). The fascinating idea of peri¬ 
odic comet impacts causing ecological catastrophies emerged 
(Alvarez & Muller 1984, Davis et al. 1984). It was suggested 
that an unseen solar companion (Nemesis) might induce grav¬ 
itational disturbances to the Oort comet cloud triggering peri¬ 
odic cometary showers (Davis et al. 1984, Whitmire & Jackson 
1984). Other astronomical models have been proposed later to 
account for the above periodicities, the “galactic carrousel” be¬ 
ing perhaps the most widely accepted model (e.g. the review 
by Rampino & Haggerty 1996, and references). The main idea 
of the “galactic carrousel” model is that the Oort comet cloud 
is periodically perturbed by galactic tides as the Solar System 
revolves around the centre of the Milky Way galaxy. 

We will show that only one extremely significant regularity 
exists in the impact crater record: the “human-signal”. 

2. Data 

We chose n = 82 impact craters with a diameter Di [km] and 
an age ti [Myr], which had an error (ntj in the database main- 
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tained by the Geological Survey of CanadaJil These data are 
only published in electronic form in our Table 1. Eight simul¬ 
taneous events were combined (^ n = 74) with the relations 
ni ,2 = min[o-ti, crt 2]/\/2 and Di ,2 = {Dl + where nti, 

at 2 , D\ and D 2 refer to the individual events. The geographi¬ 
cal coordinates of these pairs imply an occurrence of a double 
impact, except for one pair (i = 37). Combining probable dou¬ 
ble impacts leaves our result unchanged, but provides better 
statistics. Six subsamples were selected from Table 1 (n = 74): 

Ci: 5 < t 

C 2 : t < 250, (Jt < 20, D > 5 
C 3 : 5 < t < 300, at < 20 
C4: Cl, t is not a multiple of 5 
C 5 : C 2 , t is not a multiple of 5 
Ce: C 3 , t is not a multiple of 5 

Criteria similar to C 2 and C 3 have been applied earlier (Grieve 
& Pesonen 1996, Matsumoto & Kubotani 1996). We also anal¬ 
ysed the sample (C 7 : n = 13), where the 28.4 Myr cycle was 
originally detected (Alvarez & Muller 1984), and one sample 
of mass extinctions of species (Cs: n=8), recently compared to 
the impact crater record (Matsumoto & Kubotani 1996). Two 
sets of weights were derived: 

lOtj = At Di, 

where the constants = AD = n/[X;jLi A] 

ensure that X)" ^*1 ~ The two largest (uimax.i, 

Wmax, 2 ) and the smallest (wmin) weights, the ratio Wr = 
Wmax,l/Wmin and the average S = (wmax.l + Wmax,2)/2 for Cl, 
..., C7 are given in Table 2 (No at^ were available for Cs). 

3. Analysis 

Two types of techniques have been mainly applied to search 
for periodicity in the impact crater record: the power spectrum 
method and those discussed by Yabushita (1991). The 28.4 
Myr cycle was detected with the former technique (Alvarez & 
Muller 1984). Techniques of the latter type have been applied, 
e.g. by Yabushita (1991) and Grieve & Pesonen (1996). Yet 
it has not been fully realized that these techniques are most 
sensitive to uni-modal phase distributions. These phases are 
cireular data: a random sample of measurements representing 
(()i at ti folded with a period P (i.e. = FIX[ti/P], where FIX 

removes the integer part of U/P). Several nonparametric meth¬ 
ods for detecting both uni- and multi-modal (j)i distributions 

^ //gdcinfo.agg.emr.ca/cb-bin/crater/crater_table?e 
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{M=10,G=1,P„,„ = 2.20.P„„=200.00) (n = 61,7=0.001,m = 904.Zo=0.376) 



Fig. 1. The ten best period candidates detected with the K-method for Ci with n = 61: the overfilling factor (Jetsu Pelt 1996: Eq. 12) 
is [A(p]~^ = GM = 10 for this test between Pmin = 2.2 and Pmax = 200. The statistics for m = 904 independent frequencies yield the level 
P{Vn > 2 : 0 ) =7 = 0.001 outlined with a horizontal line. The numbers above the short vertical lines denote the locations of Pi,...,Pio. 


(M=10.G=l,P^i„ = 2.20,P^g^=200.00) (n = 23,/5„= 10.190,7 = 0.00 l,m= 108,C^=0.092) 


Pi = 2.50 
P2=5.00 
P3=3.00 
P,i=2.33 
P5=7.01 
P6=5.50 
P7=3.75 
Pe=14.88 
P9=2.68 
P,0=3.50 

0.00 0.04 0.08 0.12 0.16 0.20 0.24 0.28 0.32 0.36 0.40 0.44 0.48 



Fig. 2. The ten best period candidates detected with the WSD-method for C6('ir’D) with n = 23: the level P(Van ^ Cn/Z^n) = 7 = 0.001 for 
m= 108 is outlined with a horizontal line (Pmirn Pmax, G and M as in Fig. 1). 


Table 2. The weights Wt and wb of Ci, ..., C?: the two largest 
(t^inax,!, t^max,2) and the Smallest (wmin) weights, the ratio 
Wr = Wmax.i/Wmin, the average s = (wmax.i + Wmax.2)/2, and 

the breakdown parameter R{s) (Jetsu 1996: Eq. 5). 


w 

n 

Wmin 

'^max,l 

'^^max,2 

Wr 

s 

Ris) 

Ci(wt) 

61 

3.15-10““ 

50.41 

6.43 

1.60-10' 

28.42 

3.57 

Ci(wd) 61 

1 

0 

8.32 

6.93 

187.50 

7.63 

0.10 

C 2 (wt) 

34 

2.75-10“‘’ 

27.54 

4.41 

1.00-10'’ 

15.97 

4.21 

C2(wd) 34 

0.15 

4.66 

2.74 

30.91 

3.70 

0.10 

C3(wt) 

35 

TsTifr^ 

29.00 

3.70 

1.60-lO" 

16.35 

3.94 

C3(wd) 35 

6.99-10“^ 

4.75 

2.80 

68.00 

3.77 

0.10 

C4(wt) 

27 

i.ie-io-* 

22.44 

2.86 

1.94-10'’ 

12.65 

4.14 

C4(wd) 27 6.47-10“^ 

6.26 

3.55 

96.77 

4.91 

0.23 

C5(wt) 

25 

2.03'10~‘’ 

20.26 

3.24 

1.00-10'’ 

11.75 

4.30 

C5(wd) 25 

0.15 

4.66 

2.74 

30.91 

3.70 

0.17 

C6(wt) 

23 

1.20'10“'* 

19.15 

2.44 

1.60-10'’ 

10.80 

4.23 

C6(wd) 23 

8.17-10“^ 

4.48 

2.63 

54.84 

3.56 

0.19 

Criwt) 

13 

1.09'10~^ 

8.90 

1.09 

816.33 

4.99 

1.41 

Criwo) 

13 

0.15 

2.94 

2.36 

20.00 

2.65 

0.33 


exist. We applied two such methods by Kuiper (1960: the K- 
method) and Swanepoel & De Beer (1990: the SD-method), 
and their weighted versions (Jetsu & Pelt 1996: the WK- and 
the WSD-methods), which can utilize the additional informa¬ 
tion in Wt and wn. Our notations are as in Jetsu & Pelt (1996). 

The limit t > 5 eliminates a bias (Ci, C3, C4 and Ce), be- 
canse over 17 % of ti in Table 1 are below this limit, e.g. these 
ti have (f>i < 0.05 for P > 100. A limit in t is unnecessary for 
C2 and Cs, since the criteria in at and D eliminate most of the 
ti < 5. The criterion in D was applied, because earlier studies 
have indicated different periodicities for ti of larger and smaller 
craters (e.g. Yabushita 1991). 

The statistical “null hypothesis”: 


Ho'- “The 4>i ofti with an arbitrary period P are randomly 
distributed between 0 and 1.” 

was tested. All tests were performed between Pmin = 2.2 and 
Pmax = 200. The preassigned significance level for rejecting 
Ho was 7 = 0.001. Two examples of these tests are displayed 
in Figs. 1 and 2. Most of the detected P reaching 7 = 0.001 
are integers or ratios of two integers (see Table 3). We call 
this regularity arising from the bias of assigning integer val¬ 
ues for ti as the “human-signal”. Unfortunately, any arbi¬ 
trary P can be expressed by a ratio of two integers. Two 
Pi and P 2 of the “human-signal” will induce a set of spuri¬ 
ous periods P' = + ki{k 2 P 2 )~^]~^ (fci = ±1,±2,... and 

k 2 — 1, 2,...). Fnrthermore, another set of spurious periods will 
he P' — koPiP 2 (fcs = 1,2,...). The “human-signal” is strongest 
in Cl, since many larger ti end with 5 or 0, hence C 4 , C 5 and 
Ce were selected. An arbitrary Pmin = 2.2 was chosen to avoid 
the detection of, e.g. P = 2, 3/2, and 1. Note that the order of 
significance for the period candidates in Figs. 1 and 2 is not 
necessarily the same as in Table 3, since these periodograms 
were derived for the overfilling factor = GM — 10 (Jetsu 

& Pelt 1996: Eq. 12). The final values of Table 3 in the vicinity 
of these period candidates were determined with [A(/]“^ = 100. 
The four nonparametric methods are not “equally sensitive” 
to different types of </i distributions (Jetsu 1996, Jetsu & Pelt 
1996), which explains the differences in the detected P of Table 
3. 

Althongh nncertainties in some critical level estimates Qk, 
Qwk, Qsd and Qwsd exist, these do not alter the order of 
significance for the detected P (Jetsu & Pelt 1996). Thus the 
“human-signal” is the most significant periodicity in the data. 
Firstly, the estimate for the number of independent frequen¬ 
cies in all tests was m = (/max —/min)//o, where /max — ^min ’ 
/min = Pmix and /o = (tmax “ tmin) “ ^ These m estimates were 
checked with the empirical correlation function r(fc) (Jetsu 
& Pelt 1996: Eq. 14). We denote the number of indepen- 
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Table 3. The five best periods detected in Ci, Cg with the K-, WK-, SD- and WSD-methods: The estimate of the number 
of independent frequencies is m for a sample of a size n. The symbols J and f denote the cases m>m' and m<rd. The critical 
levels for each period P are Qk Qwk Qsd and Qwsd. The rejection of Ho with 7 = 0.001 is indicated by •. Finally, the cases 
when an application of some particular method would not yield reliable results are denoted by * 



K-method 
no weights 

P (Qk) 

WK-method 

Wd 

P (Qwk) 

SD-method 
no weights 

P (Qsd) 

WSD-method 

Wd 

P (Qwsd) 

Cl 

2.50 (~ Ot) • 

2.50 (~ Ot) • 

20.00 (Ot) • 

2.22..( (Ot) • 

n = 61 

5.00 (2.0'10'^®t) • 

3.00 (l.O-lO-^^t) • 

2.66... (OJ) • 

5.00 (Ot) • 

m = 904 

3.00 (3.9-10“®t) • 

5.00 (1.2-10“®t) • 

5.00 (OJ) • 

2.66... (Ot) • 


10.00 (2.1-10"'‘t) • 

15.00 (5.0-10"®t) • 

8.00 (Ot) • 

4.00 (Ot) • 


3.33 (2.1-10-^t) • 

9.97 (4.1-10-’^t) • 

2.50 (Ot) • 

20.00 (Ot) • 

C 2 

3.00 (0.0078) 

2.50 (2.6-10'“) • 

5.00 (8.9-10~"'') • 

5.00 (8.9-10~"®) • 

n = 34 

2.50 (0.11) 

7.76 (6.6-lO""') • 

2.50 (8.9-10~“) • 

2.50 (8.9-10"“) • 

m = 110 

11.67 (0.58) 

3.00 (7.4-lO’"') • 

3.00 (1.6-10“®) • 

3.00 (1.9-10“®) • 


17.50 (0.62) 

9.95 (0.0037) 

2.25 (4.0-10“®) • 

6.00 (6.7-10-®) • 


7.00 (0.73) 

12.56 (0.0046) 

6.00 (5.0-10“®) • 

4.00 (7.3-10-®) • 

C3 

2.50 (0.0023) 

2.50 (1.4-10^'") • 

5.00 (4.2-10^"“) • 

2.50 (4.2-10-"®) • 

n = 35 

3.00 (0.0053) 

3.00 (1.3-10"'‘) • 

2.66... (4.2-10~^°) • 

8.00 (4.2-10-^°) • 

m = 127 

5.00 (0.18) 

7.76 (1.5-10-^) • 

2.50 (4.2-10~i°) • 

5.00 (4.2-10-“) • 


2.86 (0.58) 

12.56 (0.0016) 

4.00 (4.2-10~^°) • 

4.00 (4.2-10-“) • 


11.62 (0.80) 

9.95 (0.0038) 

10.00 (4.2-10“^°) • 

2.66.. (4.2-10-“) • 

C4 

2.50 (0.33t) 

2.50 (2.3-10-'‘t) • 

4.00 (3.M0~'’t) • 

4.00 (3.1-10-®t) • 

n = 27 

13.20 (0.83t) 

2.98 (5.3-10-^t) • 

2.50 (3.M0-®t) • 

2.50 (3.1-10-®t) • 

m — 904 

3.00 (0.88t) 

2.70 (6.3-10-^t) • 

5.00 (3.M0-®t) • 

5.00 (3.1-10-®t) • 


2.98 (0.96t) 

7.77 (7.1'10-^t) • 

10.00 (3.M0-®t) • 

10.00 (3.1-10-®t) • 


5.50 (0.97t) 

30.00 (0.0016t) 

3.33... (3.1-10"®t) • 

3.33... (3.1-10-®t) • 

C5 

2.50 (0.17) 

7.77 (0.0011) 

2.50 (2.0-10'“) • 

2.50 (2.0-10-®) • 

n = 25 

17.43 (0.31) 

2.98 (0.0011) 

5.00 (2.0-10“®) • 

5.00 (2.0-10-®) • 

m = 110 

7.50 (0.49) 

3.71 (0.0033) 

3.00 (3.7-10“®) • 

3.00 (4.9-10-®) • 


3.00 (0.49) 

7.11 (0.017) 

3.50 (1.0-10“®) • 

3.50 (1.3-10-®) • 


14.12 (0.52) 

2.67 (0.021) 

4.00 (1.7-10“®) • 

4.00 (3.4-10-®) • 

Ce 

14.00 (0.50) 

2.98 (5.7-10^^) • 

5.00 (1.0-10^^) • 

5.00 (1.0-10-^) • 

n = 23 

7.00 (0.56) 

7.77 (0.0013) 

2.50 (1.0-10“®) • 

2.50 (1.0-10-®) • 

m = 108 

3.00 (0.61) 

3.71 (0.0018) 

3.00 (1.4-10“®) • 

3.00 (1.9-10-®) • 


3.87 (0.69) 

13.54 (0.012) 

7.00 (1.9-10“®) • 

3.50 (2.9-10-®) • 


2.50 (0.75) 

2.52 (0.013) 

3.50 (2.2-10“®) • 

3.75 (3.8-10-®) • 

C7 

20.63 (0.083) 

* 

4.00 (0) • 

* 

n = 13 

5.76 (0.11) 

* 

* 

* 

m = 87 

3.82 (0.11) 

* 

* 

* 


2.30 (0.12) 

* 

* 

* 


2.95 (0.12) 

* 

* 

* 

Cg 

2.76 (0.0031) 

* 

* 

* 

n = 8 

3.68 (0.064) 

* 

* 

* 

m = 105 

3.22 (0.13) 

* 

* 

* 


11.03 (0.31) 

* 

* 

* 


13.80 (0.39) 

* 

* 

* 


dent frequencies implied by r{k) with m'. Some tests had 
m! ^ m. The correct values for the critical level are smaller 
when m> m' (i.e. the significance is higher), while the case is 
opposite for m < m'. This problem is only present in Ci and 
C 4 containing less older than younger craters, and thus the 
inverse of tmax — tmin provides a poor /o estimate. Secondly, 
Qk > Qwk when the same period is detected with the K- and 
WK-methods (Jetsu & Pelt 1996: Eq. 30), which is due to the 
large scatter in wd (see Table 2). For example, the sum of the 
two largest wd in Ci(n = 61) is 15.25, which disrupts the statis¬ 
tics of the WK-method. However, the values of Qk are reliable 
when m — m'. Thirdly, the statistics of the WSD-method are 


quite robust even for a higher scatter of weights. But the val¬ 
ues of Qsd are Qwsd are uncertain for smaller samples, since 
the analytical estimate for the critical parameter Cn is accu¬ 
rate only for larger samples (Jetsu & Pelt 1996: Eq. 16). The 
SD-method reveals some bizarre cases (Qsd = 0= never), e.g. 
for Ct with P = 4. No random sample can contain so many 
time differences that are multiples of 4. Even if Cr contained 
a period that is apparently not due to the “human-signal”, 
such as the 28.4 Myr (Alvarez & Muller 1984), it is impossi¬ 
ble to decide whether this period is a multiple of, say 7x4. In 
fact, P = 28.79 (Qk = 0.13) is the sixth most significant period 
detected with the K-method in C 7 . 
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Table 2 shows, why no results for wt were obtained. The 
largest Wt are typically ^ n/2, and thus the statistics of the 
WK-method would be unreliable. As for the WSD-method, 
we refer to the breakdown parameter i?(s) in Table 2. If i?(s) 
exceeds unity, the statistics of the WSD-method are disrupted 
(Jetsu 1996: Eq. 5). Table 2 indicates that the WK- and WSD- 
methods must not be applied with wt ■ Because the high scatter 
in Wt prevents applications of the WK- and WSD-methods, 
we checked, if the removal of less accurate ti would significantly 
alter Wmax.i or R(s) for wt- But this removal of less accurate ti 
did not yield R(s) < 1 for Wt even if this procedure was carried 
out until only 5 values remained in Ci, ..., C 7 . The same applies 
to Wmax.i, which remains too close to n. 

Why was the “human-signal” not detected earlier, while 
the 28.4 Myr cycle was detected? The answer to the first ques¬ 
tion is that the </>i distributions connected to the “human- 
signal” are mostly multi-modal, and can not be detected with 
methods sensitive to uni-modal distributions. There are 
several answers to the second question. The 28.4 Myr cycle, 
as well as any other noninteger period, may be induced by 
the “human-signal”. For example, one earlier study revealed 
mainly multiples of 5 (Yabushita 1991: e.g. 30x1 =15x2= 10x3 
= 6x5). Since C 7 has Qsd = 0 for P = 4, no Monte Carlo simu¬ 
lation will ever “succeed” in producing so many exactly equal 
time differences, let alone the two additional equal t values in 
C 7 , i.e. the earlier significance estimates were not correct (Al¬ 
varez & Muller 1984). Finally, the scatter of wt in C 7 is so large 
that the case n = 13 does not occur, because the sum of the 
two largest Wt is 9.0. An analysis of the superimposed gaus- 
sians of these ti detects periods from the highest peaks with 
large Wt, while the smaller Wt do not influence the result (e.g. 

Wmin = 0.01). 

The small sample Cs could only be reliably analysed with 
the K-method. No signs of the 26 Myr cycle were detected 
(Raup & Sepkoski 1984, 1986), nor any period with Qk< 7 = 
0.001. The best periods do not betray any trivial signature of 
the “human-signal”, although half of the ti are integers. If the 
most significant P = 2.76 (Qk = 0.0031) represents real peri¬ 
odicity, Raup & Sepkoski (1984, 1986) have identified about 
every tenth mass extinction event. 

4. Conclusions 

A few topics must be emphasized to avoid misunderstanding, 
(i) We did not assume that the integer ti cause the detected 
periodicities. On the contrary, the periodicities were uniquely 
detected with non-parametric statistical methods (i.e. model 
independent). We simply tested the “null hypothesis” (Ho) that 
the impact crater ages represent a random sample of cireular 
data. The “human-signal” reaches 7 = 0.001 in all subsamples 
Cl, ..., C 7 , the critical levels of the SD- and WSD-methods be¬ 
ing extremely high. The analytical statistics of our methods are 
robust, i.e. Monte Carlo or other computational techniques are 
unnecessary. (ll) The “human-signal” induces irregular multi¬ 
modal (fi distributions, which were not detected earlier with 
methods sensitive to uni-modal (j>i distributions. For example, 
the power spectrum method is most sensitive to sinusoidal vari¬ 
ations. If the (pi distribution for the “correct” P' were exactly 
bi-modal, the peaks of the power spectrum would be at P', 
P' 12, ... But the power spectrum method is quite insensitive to 
more irregular (pi distributions. That the “human-signal” was 
not detected earlier, over twelve years after the study by Al¬ 


varez & Muller (1984), is a direct consequence of favouring 
methods sensitive to uni-modal (pi distributions. Why should 
the (pi distributions connected to the possible periodicity in 
terrestrial impact cratering rate or mass extinctions of species 
actually be uni-modal or of any regular shape? (ill) It may 
well be that the large ut, prevent detection of periodicity (e.g. 
Heisler & Tremaine 1989). In that case our study is simply 
an exercise of statistics providing one new argument against 
real periodicity. However, considering the prevailing theories 
based on assuming the presence of real periodicities (see e.g. 
Rampino & Haggerty 1996: “Shiva Hypothesis”), it was a high 
time to perform this exercise, (iv) The “human-signal” has 
most probably induced those spurious periods with more or 
less unimodal (pi distributions, which were detected in several 
earlier studies (e.g. Alvarez & Muller 1984, Yabushita 1991). 
In any case, the “human-signal” is clearly the most significant 
periodicity in the impact crater record, (v) We do not argue 
that the major comet or asteroid impacts and the mass ex¬ 
tinctions of species are uncorrelated, but emphasize that the 
“human-signal” dominates the time distribution of the former 
events. 

Our conclusions are simple. The epochs of mass extinc¬ 
tion events of species may follow a possibly “nonhuman” cycle 
of 2.76 Myr, but the currently available impact crater data 
definitely reveals the embarrassing “human-signal”. The fel¬ 
low scientists have unconsciously offered a helping hand to the 
Nemesis (e.g. Davis et al. 1984) or “galactic carrousel” (e.g. 
Rampino & Haggerty 1996). The arduous task for the future 
geological research is to determine more accurate (preferably 
noninteger) revised ages for impact craters to eliminate the 
“human-signal”, which may then lead to a detection of real 
periodicity. Over a decade has elapsed in redetecting the regu¬ 
larities of our own integer number system and then interpreting 
them as periodicity in the ages of impact craters. 
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Table 1: The ages {t) of the impact craters of a diameter (D) with an error estimate (crt) for t. The locations and the values of at 
and D for the combined events are given in columns 5, 7 and 8, while “Y” denotes the samples Ci, Ce- 


i Crater 


Country 


t[Myr] (Tti/^ta 73[km] D 1 /D 2 Location 


Cl C 2 C 3 C 4 C 5 Ce 


1 Kaalijarvi 

2 Wabar 

3 Boxhole 

4 Barringer 

5 Lonar 

6 Pretoria Saltpan 

7 Zhamanshin 

8 Bosumtwi 

9 New Quebec 

10 Tenoumer 

11 Aouelloul 
ie~El’gygytgyn 
l^&oter Kamm 
laffiigach 
l ygh unak 
levies / Steinheim 
I'^Haughton 
ISicgancha 
l)i?opigai 
20Chesapeake Bay 
2lOVanapitei 
234^istastin 

23 Logoisk 

K aiyli 

amensk/Gusev 
ontagnais 
agozinka 
arquez 

2#^hicxulub 
3ft^ara/Ust-Kara 
'lanson 
jappaiarvi 
5oltysh 
Dellen 
iteen River 
Deep Bay/West Hawk 
tijswellV Zapadnaya 
?eleny Gai 
Mien 

_ fcokoonooka 
4 Ic^omistrovka 
42:Cosses Bluff 
43 Mjplnir 
44i?:jiverpool 
4S^epriai 

4)^^uchezh-Katunki 
4'4-Rochechouart 
4gt^ells Creek/Red Wing 
49 Manicouagan 
50 Obolon’ 

51 Saint Martin 
52 Araguainha Dome 
53 Kursk 

54 Clearwater West/East 

55 Dobele 

56 Crooked Creek 

57 Charlevoix 

58 Flynn Creek 

59 Siljan 

60 Kaluga 

61 Ilyinets 

62 La Moinerie 

63 Couture 

64 Calvin/Brent 

65 Pilot 

66 Ames 

67 Gardnos 

68 Holleford 

69 Mizarai 

70 Janisjarvi 
71 Morokweng 
72 Teague 
73 Sudbury 
74 Vredefort 


Estonia 
Saudi Arabia 
Australia 

U.S.A. 

India 

South Africa 

Kazakhstan 

Ghana 

Canada 

Mauritania 

Mauritania 

Russia 

Namibia 

Kazakhstan 

Kazakhstan 

Germany 

Ganada 

Russia 

Russia 

U.S.A. 

Ganada 

Ganada 

Belarus 

Kazakhstan 

Russia 

Ganada 

Russia 

U.S.A. 

Mexico 

Russia 

U.S.A. 

Finland 

Ukraine 

Sweden 

Ganada 

Canada 

Canada/Ukraine 

Ukraine 

Sweden 

Australia 

Ukraine 

Australia 

Norway 

Australia 

Lithuania 

Russia 

France 

U.S.A. 

Canada 

Ukraine 

Canada 

Brazil 

Russia 

Canada 

Latvia 

U.S.A. 

Canada 

U.S.A. 

Sweden 

Russia 

Ukraine 

Canada 

Canada 

U.S.A./Canada 

Canada 

U.S.A 

Norway 

Canada. 

Lithuania 

Russia 

South Africa 

Australia 

Canada 

South Africa 


0.004±0.001 

0.006±0.002 

0.03±0.0005 

0.049±0.003 

0.052±0.006 

0.22±0.052 


0.11 

0.097 

0.17 

1.186 

1.83 

1.13 


0.9±0.1 


13.5 



Y 



Y 


1.03±0.02 


10 5 



Y 



Y 


1.4±0.1 


3.44 








2.5±0.5 


1.9 








3.1±0.3 


0.39 








3.5±0.5 


18 



Y 



Y 


3.7±0.3 


2.5 








6±3 


7 


Y 

Y 

Y 

Y 

Y 

Y 

12±5 


3.1 


Y 


Y 

Y 


Y 

15±0.7 

(1/1) 

24.3 

(24/3.8) N48.9E10.6/N48.0E10.1 

Y 

Y 

Y 




23±1 

24 

Y 

Y 

Y 

Y 

Y 

Y 

25±20 


20 


Y 

Y 

Y 




35±5 


100 


Y 

Y 

Y 




35.5±0.6 


85 


Y 

Y 

Y 

Y 

Y 

Y 

37±2 


7.5 


Y 

Y 

Y 

Y 

Y 

Y 

38±4 


28 


Y 

Y 

Y 

Y 

Y 

Y 

40±5 


17 


Y 

Y 

Y 




46±7 


5.5 


Y 

Y 

Y 

Y 

Y 

Y 

49±0.14 

(0.2/0.2) 

25.2 

(25/3.5) N48.3E40.2/N48.4E40.2 

Y 

Y 

Y 

Y 

Y 

Y 

50.5±0.76 

45 

Y 

Y 

Y 

Y 

Y 

Y 

55±5 


9 


Y 

Y 

Y 




58±2 


13 


Y 

Y 

Y 

Y 

Y 

Y 

64.98±0.05 


170 


Y 

Y 

Y 

Y 

Y 

Y 

73±2.1 

(3/3) 

69.6 

(65/25) N69.2E65.0/N69.3E65.3 

Y 

Y 

Y 

Y 

Y 

Y 

73.8±0.3 

35 

Y 

Y 

Y 

Y 

Y 

Y 

77.3±0.4 


23 


Y 

Y 

Y 

Y 

Y 

Y 

88±3 


24 


Y 

Y 

Y 

Y 

Y 

Y 

89±2.7 


19 


Y 

Y 

Y 

Y 

Y 

Y 

95±7 


25 


Y 

Y 

Y 




100±35 

(50/50) 

13.2 

(13/2.44) N56.4W103.0/N49.8W95.2 Y 






115±7.1 

(10/10) 

39.2 

(39/4) N58.4W109.5/N49.7E29.0 

Y 

Y 

Y 




120±20 

2.5 

Y 


Y 




121±2.3 


9 


Y 

Y 

Y 

Y 

Y 

Y 

128±5 


55 


Y 

Y 

Y 

Y 

Y 

Y 

140±20 


2.7 


Y 


Y 




142.5±0.5 


22 


Y 

Y 

Y 

Y 

Y 

Y 

143±20 


40 


Y 

Y 

Y 

Y 

Y 

Y 

150±70 


1.6 


Y 






160±30 


8 


Y 






175±3 


80 


Y 

Y 

Y 




186±8 


23 


Y 

Y 

Y 

Y 

Y 

Y 

200±18 

(100/25) 

15 

(12/9) N36.4W87.7/N47.6W103.6 Y 

Y 

Y 




214±1 

100 


Y 

Y 

Y 

Y 

Y 

Y 

215±25 


15 


Y 






220±32 


40 


Y 






247±5.5 


40 


Y 

Y 

Y 

Y 

Y 

Y 

250±80 


5.5 


Y 






290±14 

(20/20) 

44.4 

(36/26) N56.2W74.5/N56.1W74.1 

Y 


Y 




300±35 

4.5 

Y 






320±80 


7 


Y 






357±15 


54 


Y 



Y 



360±20 


3.55 


Y 






368±1.1 


52 


Y 



Y 



380±10 


15 


Y 






395±5 


4.5 


Y 






400±50 


8 


Y 






430±25 


8 


Y 






450±7.1 

(10/30) 

9.3 

(8.5/3.8) N41.8W86.0/N46.1W78.5 

Y 






455 ±2 

6 

Y 






470±30 


16 


Y 






500±10 


5 


Y 






550±100 


2..35 


Y 






570±50 


5 


Y 






698±22 


14 


Y 



Y 



1400±200 


70 


Y 






1630±5 


30 


Y 






1850±3 


250 


Y 






2018±14 


300 


Y 



Y 







